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Abstract 

We determine families of spherically symmetrical /^-dimensional quantum potential functions 
V(r) having ground state wavefunctions that exhibit, either in configuration or in momentum space, 
the form of an isotropic g-Gaussian. These wavefunctions admit a maximum entropy description in 
terms of S q power-law entropies. We show that the potentials with a ground state of the g-Gaussian 
form in momentum space admit the Coulomb potential — 1/r as a particular instance. Furthermore, 
all these potentials behave asymptotically as the Coulomb potential for large r for all values of the 
parameter q such that < q < 1. 
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I. INTRODUCTION 



Extended versions of the maximum entropy principle based upon power law S q entropies 
[TJ [2] have been found to provide useful tools for the description of several physical systems 
or processes [THE]. Indeed, various important equations in mathematical physics admit 
exact solutions of the maximum S q form such as, for example, the polytropic solutions 
to the Vlasov-Poisson equations j3], time dependent solutions to some evolution equations 
involving non linear, power-law diffusion terms [HE], or stationary phase-space distributions 
for Liouville equations describing anomalous thermostatting processes [6j. 

The application of information theoretical ideas to the study of the eigenstates of diverse 
quantum systems has attracted the attention of researchers in recent years [7HH] . The stan- 
dard maximum entropy principle, based on the optimization of Shannon's entropic measure 
under appropriate constraints, plays a distinguished role within these lines of enquiry. This 
principle has been successfully applied to the characterization of the eigenstates of various 
quantum systems (see, for instance, [T3HT5] and references therein). In particular, it is well 
known that the probability densities in both position and momentum space corresponding to 
the ground state of the isotropic D-dimensional quantum harmonic oscillator are Gaussians, 
which are probability densities maximizing Shannon' entropy under the constraints imposed 
by normalization and the expectation value of the square r 2 of the radial coordinate. 

It would be of considerable interest to extend to the S^-based framework the maximum 
entropy approach to the description of the eigenstates of quantum systems. This formalism 
has already been applied to the study of various quantum phenomena (see, for example 
[T61IPT] ). but its application to characterize the probability densities associated with quantum 
eigenstates remains largely unexplored. The maximum entropy formalism based on the S q 
entropies leads to a generalization of the Gaussian probability density, which is given by 
the so-called g-Gaussians [TJ [2]. These g-Gaussian constitute some of the simplest and 
most important examples of maximum-S , q distributions. The aim of the present work is to 
determine the form of those spherically symmetric quantum potentials V(r) whose ground 
state wavefunctions (in position or in momentum space) are associated with g-Gaussian 
densities. 
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II. Q-GAUSSIAN GROUND STATES IN CONFIGURATION SPACE 



We are going to consider a spinless particle of mass m in a D-dimensional configuration 
space. The eigenfunctions ip(r) associated with a potential V(r) obey then the Schrodinger 
equation, 

- |^V 2 ^ + ViP = (1) 

where V 2 is the D-dimensional Laplacian operator, h is Planck's constant and E is the 
energy eigenvalue. We assume in the rest of this paper that m = h = 1. Since we are going 
to consider spherically symmetric potentials, the Schrodinger equation for the concomitant 
ground states (which are spherically symmetric) simplifies to 

1 9 f r o-m +vi , =E ^ (2) 



2 r D-l Q r ^ Q r 

where 

' D \ 1/2 

£*? ( 3 ) 

<i=l J 

is the radial coordinate. 

Let us consider a D— dimensional spherical q— Gaussian wavefunction in the configuration 
space 

V>(r) = C (1 - (q - 1) f3r 2 )^ (4) 

where q and /3 are positive parameters and C is an appropriate normalization constant. If 
q < 1 the g-Gaussian wavefunction Q remains finite for all r e M 15 . On the other hand, when 
g > 1 the g-Gaussian vanishes at r = 1/ \/{q — and is set to zero for r > 1/ yjjq — l)/3 
(see below a discussion on the physical meaning of this cut-off). The space probability 
density p(r) = \ip (r) | 2 associated with the wavefunction Q maximizes Tsallis' power-law 
entropic functional 

under the constraints given by normalization and the expectation value of r 2 |;2J (it can 
also be regarded as a probability density maximizing Renyi's functional under the same 
constraints). 



Replacing Q into (JTj) and after some algebra, we find that the wavefunction is the ground 
state of the potential 



V 



B \-D + Br 2 (D (q - 1) + 3 - 2q) 
2 [ (1 - (q- I) (Br 2 f 



(6) 



with eigenenergy equal to 0. 

In the case q < 1 the potential function ([6j is finite for all r G 1R D . On the other hand, 
when q > 1 the potential function is singular when r adopts the particular value 



Physically, this means that when q > 1 the potential function ([6]) has an "infinite wall" at 
r = r w and the quantum particle is confined within the region r < r w . In this case, the 
g-Gaussian wavefunction ^ vanishes at r = r w and must be set equal to zero when r >r w . 
This constitutes an example of the so-called Tsallis' cut-of condition [21 E]- 

In the limit q — v 1 the g-Gaussian wavefunction Q becomes a standard Gaussian, and the 
potential function ^ reduces to the D-dimensioanl isotropic harmonic oscillator potential 
(notice that the origin of the energy scale is shifted) 



The one dimensional instance of the potential (|6| has been studied in |18j . This potential 
exhibits the interesting feature of approximate shape invariance (see [18] for details). This 
approximate symmetry becomes exact in the limit q — y 1. 

III. Q-GAUSSIAN GROUND STATES IN MOMENTUM SPACE 

We now look for solutions solution of the Schrodinger equation having the form of a 
q— Gaussian in momentum space 




(7) 



V{r) 




(8) 



4> (p) = C (1 - (q - 1) (Bp 2 ) 2 <'- 1 ' 



(9) 



where 
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p 2 = X> 2 . (io) 



1=1 



As in the previous case of g-Gaussians in configuration space, q and are positive parameters 
and C is a normalization constant. We are going to consider g-Gaussians in momentum space 
with q < 1. Our aim is to determine potential functions V(r) having a ground state that, 
in momentum space, has the form (|9]). In order to do this, it will prove convenient not to 
work directly with the Schrodinger equation in momentum space but, instead, to determine 
first the Fourier transform if) (r) of if> (p) and then to consider Schrodinger's equation in 
configuration space. 

The Fourier transform of the q— Gaussian wave function ^ is 

^ (r) = wf Kv (r) 

where K v is the Bessel function of the second kind and r = |r|. The parameter parameter v 
is given by 

—l-W^Ty < 12 > 

Theorem. The function ip u (r) is solution of the Schrodinger equation associated with a 
potential 

v ^ = -2{ 1 + w^T))^W (13) 

As a special case, when the parameter v = d + \ is half integer, this potential is of the 
form 

d+ * K > 2 V 2d - 1 J Pd (r) K 1 

where p d (r) is the Bessel polynomial of degree d. 

Proof. The derivation rule for the function if} (r) is 

l -^ u (r) = -—^—^ 1 {r) (15) 
r or 2 [v — 1) 

so that the Laplace operator reads 



5 



The first term is 



1 d ( _j3_ x dip u 
dr 



r 



dr 



dr 2 



(16) 



and the second term 



d 2 ^ 



1 d 



dr 2 2 {y — 1) dr 



2(i/-l) 



(17) 



1 



-VV-i (r) + 



1 



2(z/-l) r ^ lv y 4 (i/ — 1) (i/ - 2) 



rV,- 2 (r) (18) 



so that the Laplace operator applied to ip v is 



1 d ( u^&lpv 



D 



(r) + 



1 



r D-i dr y dr j y 2(z/-l) r v ' 4 (i/ - 1) (i/ - 2) 
Moreover, the function Bessel function if w obeys the difference equation 



r a ^-a (r) (19) 



so that 



and 



We deduce 



tK v (r) = rK u ^ 2 (r) + 2 (z/ - 1) K u ^ (r) 



r"^ (r) = r z r v - 2 K v _ 2 (r) + 2 (i/ - 1) r v ~ Y K v _ x (r 



(r) = r 2 



1 



4(i/-l) (u-2) 



ipu-2 (r) + (r) 



AV>„ (r) 



1 5 f r D-i9^ 



D 



2(i/-l) 



Consequently, 



(20) 



(21) 



(22) 



1 VV-i (r) + ^ (r) . (23) 



- -A^, (r) + 



I / '., D \ V,-i(r) 



2 V 2{u-l)J ^(r) 
which means that ip u (r) is an eigenfunction of the potential V u (r) given by equation (13), 

with eigenvalue equal to — |. 

□ 
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IV. SPECIAL CASES AND ASYMPTOTICS 



A. Asymptotics 



The asymptotics for large r of the potential (13) can be computed using [T9l 9.7.2] 



*.<r)~J^(l + ^ + ...] (25) 



so that the asymptotics for the potential (13) reads 



We see then that, for large values of r, the asymptotic behavior of the potential V u (r) is 
dominated by a Coulomb-like term. 



B. Special cases 

1. Coulomb potential: taking v = | and remarking that if)i (r) = exp (— r) and t/>_i (r) = 
— fyl ( r ) , we deduce that 

" (r) - (r) = (r) (27) 

which is the Schrodinger equation associated with a Coulomb potential. The associated 
probability density in configurational space can be obtained as the squared modulus 
of the inverse Fourier transform of the ground state wavefunction in momentum space, 

^(p)oc(l + |p| 2 r^. (28) 

The momentum space representation of the eigenfunctions corresponding to the D- 
dimensional Coulomb potential have been studied in detail by Aquilanti, Cavalli and 
Coletti in [20]. 

The q- value characterizing the ground state of the — - potential is different from one. 
Indeed, it depends on the value of the space dimension D, 

« = DTP < 29 > 
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l+D 
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"(3+D)(r+l)" 
r 2 +3r+3 


1 
2 


(5+.D)(r 2 +3r+3) 
r 3 +6r 2 +15r+15 



TABLE I: Values of the potential for different half-integer values of the parameter v 



2 4 6 8 10 12 14 16 18 20 





FIG. 1: The potential functions V v appearing in Table I, corresponding to v equal to 
(top to bottom), for D = 1 (solid line) and D = 3 (dashed line). 



|, | and \ 



2. In Table I we give a few potentials resulting from different half integer values of v. 
These potentials are depicted in Fig.l. If v — -= + d, with d integer, then the entropic 
parameter q characterizing the g-Gaussian is given by 

D + 2d 



Q 



D + 2d+l 



(30) 



It is interesting that, for a given fixed value of d, we have that q — > 1 when D — > oo. 
That is, when the space dimension tends to infinity the q- Gaussian describing the 
ground state in momentum space approaches a standard Gaussian. 



V. CONCLUSIONS 

We have determined the D-dimensional spherically symmetric potential functions having 
ground states of the g-gaussian form, either in configuration or in momentum space. In the 
case of g-gaussian ground states in configuration space we obtained a bi-parametric family 
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of potentials admitting the D-dimensional isotropic harmonic oscillator as the particular 
case corresponding to the limit q — > 1. On the other hand, when considering ground states 
having the shape of a q-Gaussian in momentum space we obtained a family of potentials 
closely related to the D-dimensional Coulomb (or Hydrogen) potential — -. In point of 
fact, this family admits the standard (D-dimensional) Coulomb potential itself as a partic- 
ular instance. Moreover, for large values of r, all the above mentioned potentials behave 
asymptotically as — £ for all < q < 1. 

Within classical mechanics, it is already well known that there is a close relationship 
between the potential function — - (describing Newtonian gravitation) and maximum S q 
distributions. The celebrated polytropic solutions of the Vlasov-Poisson equations, widely 
used in the study of self-gravitating astrophysical systems, have indeed the S^-maxent form, 
and the associated velocity distributions are g-Gaussians. It is an intriguing fact that, as 
we have shown in the present work, there also exits a close connection between g-Gaussians 
and the — £ potential in quantum mechanics. 
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